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ABSTRACT 

Recent simulations of cluster formation have demonstrated that condensation of baryons into central galaxies 
during cluster formation can drive the shape of the gas distribution in galaxy clusters significantly rounder out 
to their virial radius. These simulations generally predict stellar fractions within cluster virial radii that are 
^ 2-3 times larger than the stellar masses deduced from observations. In this paper we compare ellipticity 
profiles of simulated clusters performed with varying input physics (radiative cooling, star formation, and 
supernova feedback), to the cluster ellipticity profiles derived from Chandra and ROSAT observations, in an 
effort to constrain the fraction of gas that cools and condenses into the central galaxies within clusters. We 
find that local relaxed clusters have an average ellipticity of e = 0.18 ±0.05 in the radial range of 0.04 < 
r / r 500 < L At larger radii r > O.lrsoo, the observed ellipticity profiles agree well with the predictions of 
non-radiative simulations. In contrast, the ellipticity profiles of simulated clusters that include dissipative gas 
physics deviate significantly from the observed ellipticity profiles at all radii. The dissipative simulations 
overpredict (underpredict) ellipticity in the inner (outer) regions of galaxy clusters. By comparing simulations 
with and without dissipative gas physics, we show that gas cooling causes the gas distribution to be more 
oblate in the central regions, but makes the outer gas distribution more spherical. We find that late-time gas 
cooling and star formation are responsible for the significantly oblate gas distributions in cluster cores, but 
the gas shapes outside of cluster cores are set primarily by baryon dissipation at high redshift (z > 2). Our 
results indicate that the shapes of X-ray emitting gas in galaxy clusters, especially at large radii, can be used to 
place constraints on cluster gas physics, making it potential probes of the history of baryonic cooling in galaxy 
clusters. 

Subject headings: galaxies: clusters: general - X-rays: galaxies: clusters 



1. INTRODUCTION 

Dark matter halos formed within the prevailing cold dark 
matter (CDM) par adigm are expected to collapse from triax- 
ial density peaks (Doroshkevich 1970; Bar deen et al.lll986l) 
and thus are generally expected to be triaxial (e.g jFrenk et al. 
119881: iDubinski & Carlberglll99U IWarren et alj[l992h . Halo 
shapes can be approximated by an ellipsoid uniquely de- 
termined by three principal axes, with lengths and orienta- 
tions that may vary as a function of distance to the halo 
center, as was dem onstrated in detail by a number of nu- 
merical simulation s dCole & Lacevll996l:lThomas et al J 1998 b 
Jing & Sutol 120021: ISuwa et al.1 12003] iHopkins et all 1 2005^ 
Kasun & Evrardl 120051: lAllgood et al-lhoOQ: iBett et al l 20oi 
Gottlober & Yepesl 120071: iPaz et al.1 (20081 Lau et al.1 Wil 
Schneider et alj 1201 lb . The mean axis ratios for halos of 
a given mass are sensitive to cosmological parame ters, par- 
ticularly the normalizatio n of the power spectrum dHo et alj 
l2006t lAllgood et aT] [2006). Observational confirmation of the 
halo triaxiality predicted by simulations therefore serves as 
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an important test of the CDM paradigm. Furthermore, the 
shape of the intracluster medium (ICM) can also affect mass- 
observable scaling relatio ns on the order of ~ 10% (e.g., 
iBuote & Humphrey! 120121) . Hence understanding and con- 
straining the triaxiality of the gas distribution is important for 
understanding sources of scatter in cluster mass-observable 
relations. 

Gravitational lensing is a potentially powerful method to 
directly probe t he shapes of halos. Recent analyses based on 
weak lensing dEvans & Bridlell2009|: lOguri et al.ll20ToL | 2012[) 
and strong lensing (Ric hard et al.l 2010; Q guri et al.1 120121) 
suggest that the shape of the total mass distribution is gen- 
erally consistent with the predictions of CDM simulations. 
However, the lensing constraints on cluster shapes are often 
confined to the central regions of clusters (in the case of strong 
lensing), or are based on the assumption that the ellipticity is 
constant with radius inside the halo (in the case of weak lens- 
ing). Stacking procedures to estimate halo shapes from shear 
signals require prior knowledge of the orientation of the axes 
of the halos, which introduces additional uncertainties in the 
estimated halo shapes when one assumes alignment between 
the observed g alaxy distribution and the underlying dark mat- 
ter dUISoIl). 

An alternative way to study halo shape is to mea- 
sure the shape of the hot X-ray emitting gas in clus- 
ters under the assumption that gas is in hydrostatic equi- 
librium an d thus follows the equ i potential surfaces of 
the cluster dBinney & Strimpell 119781; iFabricant et al.l Il984t 
IBuote & Tsaill 19951) . Given the relatively higher spatial reso- 
lution of X-ray maps compared to the lensing shear maps, and 
the fact that the measured ellipticity is independent of mea- 
sured cluster mass and concentration, we can use gas shape 
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to constrain the cluster ellipticity with higher accuracy. This 
can allow us to constrain the radial depen dence of ellipticity 
predicted by numer ical simulations (e.g. iJing & Sutol 120021: 
Allgood et al. 2006). Recent work suggests that the observed 
X-ray shap e is consistent with the prediction of the CDM 
model (e.g. jFlores et alfeOOTt iKaw ahara 2010|). 

Furthermore, the shape of dark matter distribution is quite 
sensitive to the degree of central concentration of mass. 
As baryons condense toward the center to form a cen- 
tral galaxy within a h alo, the dark matter distribution be- 
comes more spheri cal) Katz & Gu rinl Il991t Dubinski | 1994t 
Evrar detaLl 119941: iTiss era & Dominguez-Ten reirol 1 1 998b 
Kazantzidis et alj|2004t lHavashi et alj|2007t iDebattista et all 
2008; Tiss era et al.l l2010). This effect increases with decreas- 
ing radius, but is substa ntial even at half of the virial radius 
dKazantzidis et al.l 120041) . The main mechanism behind this 
effect lies in adiabatic changes of the shapes of particle orbits 
in response to the more centra lly-concentrated mass distribu- 
tion a fter baryon dissipation (Dubinski 1994; Debattist a et alj 
2008). The change in the overall shape of the dark mat- 
ter distribution m akes the shape of the ICM rounder as well 
dLau et al.H201 11) . Therefore, studies of the shape profile can 
constrain the amount of cooling that has occurred during clus- 
ter f ormation. 

In lLau et"ail d201 ll hereafter Paper I), we used a set of 16 
high-resolution hydrodynamical cosmological cluster simula- 
tions with and without dissipative gas physics and showed that 
the shape and normalization of the X-ray ellipticity profile 
are sufficiently sensitive to distinguish between the shape of 
ICM predicted in simulations with and without cooling with 
high significance. In particular, we showed that the baryon 
condensation due to gas cooling makes gas outside of clus- 
ter cores (0.1 < r/r^Q < 1), Q more spherical compared to 
non-radiative simulations. At small radii (r/rsoo < 0.1), gas 
cooling leads to a more oblate gas shape due to the formation 
of a rotating gaseous slab. In the non-dissipative simulations, 
on the other hand, the gas shape becomes more spherical in 
the central regions due to the enhanced random gas motions 
in the region with a flat entropy profile. Thus, by measuring 
the gas shape in the core, one can infer the level of gas mo- 
tions in cluster cores. By analyzing mock Chandra images of 
the simulated clusters, we showed that the baryonic effects on 
the ellipticity profiles of galaxy clusters should be discernible 
using X-ray observations. 

In this companion paper we compare a sample of local clus- 
ters observed by Chandra and ROSAT to the same suite of hy- 
drodynamical cluster simulations used in Paper I to constrain 
the fraction of gas that cooled to form the central galaxy in 
real cluster s. A similar comp arison has recently been car- 
ried out by IFang et al.1 ((2009). Our study is different from 
iFang et aTJ (120091) in two important respects. First, we an- 
alyze and compare significantly larger samples of both ob- 
served and simulated clusters. Second, we focus on the effects 
of baryon condensation on the elliptici t y of e quipotential sur- 
faces at larger radii, while IFang et aTl (|2009) focused on the 
effects of more recent cooling in the central regions of clus- 
ters. 

This paper is organized as follows. We describe the simu- 
lations in Section |2] and the X-ray data in Section [3] In Sec- 
tion|U we present the comparison between the ellipticity pro- 
files of observed and simulated clusters. We conclude with a 

7 Here f-500 is the radius of the halo where the average density within this 
radius is 500 times the critical density of the universe. 



summary and discussion in Section[5] 

2. SIMULATIONS 
2.1. Simulation data 

We use high-resolution cosmological simulations of 16 
cluster-sized systems in the flat concordance ACDM model: 
Cl m = 1 -fi A = 0.3, fl b = 0.04286, h = 0.7 and cr 8 = 0.9, where 
the Hubble constant is defined as lOO/i km s" 1 Mpc -1 , and a% 
is the power spectrum normalization on an 8/i _I Mpc scale. 
The simulations were performed with the A daptive Refine- 
ment Tree (ART) jV-b ody+gasdynamics code (iKravtsovl 19991: 
iRravtsov et al.ll2002h . an Eulerian code that uses adaptive re- 
finem ent in space and tim e, and (non-adaptive) refinement in 
mass (Klypinetal. 2001) to reach the high dynamic range 
required to resolve cores of halos formed in self-consistent 
cosmological simulations. In order to assess the effects of 
gas cooling and star formation on the cluster shapes, we con- 
ducted each cluster simulation with two different prescrip- 
tions for gas dynamics. In one set of runs we treated only the 
standard gas dynamics for the baryonic component without 
either radiative cooling or star formation. We refer to these 
as non-radiative (NR) runs. In the second set of runs, we in- 
cluded gas, metallicity-dependent radiative cooling, star for- 
mation, and superno va feedback ( CSF). Th e simulations are 
discussed in detail in Na gai et aTl (l2007albl) and we refer the 
reader to these papers for more details. Here we summarize 
the main parameters of the simulations. 

Simulations were run using a uniform 128 3 grid and 
eight levels of mesh refinement in computational boxes of 
120/T 1 Mpc and 80/T 1 Mpc on a side, respectively. These 
simulations achieve a dynamic range of 32,768 and a for- 
mal peak resolution of w 3.66/T 1 kpc and 2.44 h kpc, corre- 
sponding to an actual resolution of « 7 h~ l kpc and 5h~ l kpc 
for the 120 h~ l Mpc and 80/T 1 Mpc boxes, respectively. Only 
regions of ~ 3- 10/T 1 Mpc around each cluster was adap- 
tively refined, while the rest of the volume was followed on 
the uniform 128 3 grid. The mass resolution, m part , corre- 
sponds to an effective 512 3 particles in the entire box, or a 
Nyquist wavelength of An y = 0.469 /i -1 and 0.312/i _1 comov- 
ing Mpc for box sizes of 120 h~ l Mpc and 80/i _1 Mpc, respec- 
tively. These correspond to 0.018/i _1 and 0.006/i _1 Mpc in 
physical units at the initial redshifts of the simulations. The 
dark matter (DM) particle mass in the regions around each 
cluster was 9.1 x lO 8 /z _1 M and2.7 x 10 8 /? _1 M Q for the two 
box sizes, while other regions were simulated with lower mass 
resolution. 

2.2. Mock Chandra images 

To compare our simulation results to observations of X-ray 
clusters, we analyze mock Chandra images for each simu- 
lated cluster, with the instrumental response of real X-ray in- 
struments taken into account, and measured the ellipticity pro- 
files of the gas emission in these mock images. Below we give 
an overview of the methods used to generate the mock images. 
Detailed descri ptions of these mo ck images can be found in 
Section 3.1 of iNagai et all (l2007al) . 

First we created X-ray flux maps of the simulated clusters 
viewed along three orthogonal projections. A flux map is 
computed by projecting the X-ray emission of hydrodynamic 
cells enclosed within three virial radii of a cluster along the 
line of sight. The X-ray emissivity in each computational grid 
cell is computed as a function of proton and electron densi- 
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ties, gas temperature and metal abundance. Emission from 
gas with temperature less than 10 5 K is excluded as it is be- 
low the Chandra bandpass. We then convolved the emission 
spectrum with the response of the Chandra front-illuminated 
CCDs and draw a number of photons at each position and 
spectral channel from the corresponding Poisson distribution. 
Each map has an exposure time of 100 ks (typical for deep 
observations) and includes a background with the intensity 
corresponding t o the quiescent backgro und level in the ACIS- 

1 observations (Mar kevitch et alj|2003l) . The resolution of all 
the maps is 6 kpc pixel" 1 . 

From these data, we generated X-ray images in the 0.7- 

2 keV band and used them to identify and mask out all de- 
tectable small-scale clumps, as is routinely done in observa- 
tional analyses. Our clum p detection is based on a wavelet de- 
composition algorithm ( Vikhl inin et al.1 11998). The holes left 
by masking out substructures in the photon map are filled in 
by the values from the decomposed map of the largest scale in 
wavelet analysis. We have tested that this method preserves 
the global shape of the photon distribution well. The back- 
ground is removed when estimating ellipticities as it can bias 
them low at radii where background dominates the intrinsic 
emission. Throughout this paper we assume that the cluster 
redshift is z b s = 0.06 for the z = sample. 

To investigate how the dynamical states of clusters affects 
their shapes, we divided our sample into relaxed and unre- 
laxed sub-samples based on visual examinations of their mock 
X-ray images. Relaxed clusters are characterized by regular 
shapes and the absence of prominent substructures in all three 
of their orthogo nal mock maps. De tails of the classification 
can be found in Na gai et afl (l2007al) . We focused our com- 
parison with observations on relaxed clusters in this paper. 
There are a total of 2 1 mock images in which the clusters are 
relaxed. 

We esti mated axis ratios follow ing iDubinski & Carlbergl 
(119911) and Kazantzi dis et alJ ((2004) where we calculated the 
inertia tensors of the intracluster gas distribution within ra- 
dial shells. This method provides a robust measure of shape 
(see detailed discussion and tests by Zemp et al.ll20llb . The 
ellipticities are measured by computing the inertia tensor of 
the photon distribution weighted by the photon counts and are 
estimated within radial annular bins. Details about the shape 
estimation method can be found in Sections 3.1 and 4.1 of 
Paper I and Section l3~2l of the current paper. 

3. OBSERVATIONS 
3.1. Data 

We analyze Chandra and ROSAT/PSPC X-ray observations 
of nearby galaxy clusters. The X-ray maps used are flat- 
fielded and with point source excluded. The c luster sample 
is taken from Table 2 of Vikhl inin et alJ (l2009al) . and we refer 
the reader to this paper for a detailed discussion of selection 
and data reduction for this sample. Information on the clus- 
ter sample used in the current paper is also presented in Ta- 
ble Q] The radius rspp is der i ved fro m the best fit M500 - Yx re- 
lation from Vikhlini n et alJ (|2009a), where we assume a con- 
cordance ACDM model with f2 m = 0.3, ft\ = 0.7, and Hubble 
parameter h = 0.72. Clusters that exist in the Chandra or the 
ROSAT samples are labeled 'C or 'R' respectively in the sec- 
ond column in Table [TJ 

We use both Chandra and ROSAT data whenever possible 
to probe different radial ranges within clusters. Chandra mea- 
surements provide exquisite cluster shape measurements in 
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FIG. 1 . — Distribution of the shortest distance from the cluster center to the 
chip edge of our sample of Chandra clusters. 



the inner regions of clusters. While having a very good spa- 
tial resolution of ~ 4" (after binning to achieve good signal- 
to-noise), the ellipticity measurements are limited to within 
small radial range (r <C ^500) due to the chip geometry of the 
ACIS camera on Chandra. The centers of most of the Chan- 
dra clusters lie close to the edge of the ACIS-I chip, thus limit- 
ing the radial ranges available for robust shape measurement. 
Figure Q] shows the distribution of the shortest distance from 
the cluster center to the closest edge of the ACIS-I chip, where 
we define the edge of the chip to be the pixels that have less 
than 40% of the maximum value in the exposure maps. We 
further restrict our ellipticity measurements for each of the 
Chandra clusters to within 90% of this distance. 

The large field of view (FOV) of ROSAT enables us to ex- 
tend the ellipticity measurement to larger radii, although its 
spatial resolution is considerably worse than that of Chandra. 
The point-spread-function of ROSAT/PSPC has an FWHM of 
~ 25". Hence we limit our ROSAT ellipticity measurement 
with the innermost radial bin being greater than this value and 
only uses pixels which have fluxes above the background flux 
determined as the mean flux outside 7*500 of the cluster. 

We also check our cluster sample for any bubbles or jet-like 
structures due to active galactic nuclei (AGNs) in the core re- 
gion that could bias our ellipticity measurements. In general, 
these structures have little effect on the ellipticity, neverthe- 
less we exclude the regions that contain these structures from 
our analyses wherever they are visually present in the images. 

3.2. Methods 

We estimate the ellipticity of the cluster in an X-ray map by 
computing the inertia tensor of the X-ray photon distribution 
in annular bins. The inertia tensor is computed as 



= V 



(D 



using pixels with photon count w a if the ellipsoidal distance 
of the pixel 



(2) 
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TABLE 1 

Properties of Observed Sample of z < 0.1 Clusters 
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is within a spherical annulus [R,R+ A/?], where R is the radius 
of the annulus . The primed coordinates are the coordinates 
of the pixel rotated to the frame of the principal axes of the 
inertia tensor. The ellipticity is initially set to zero, and then 
computed iteratively as one minus the ratio of the square root 
of eigenvalues from diagonalizing the inertia tensor. The el- 
lipticity, e, is defined as 

e=l-~, (3) 
a 

where a and b are the semi-major and the semi-minor axis 
lengths from the fitted ellipse, respectively. The iteration stops 
when the ellipticity converges to within 1% of its previous 
value. 

Following Paper I, we identi fied point sources and substruc- 
tures using wavelet algorithm ( Vikhlinin et al.|[l998r) . 



As discussed in Section 2, the substructures are identified 
using wavelet algorithm (Vikhlinin et al 1998). In relaxed 
clusters only a few substructures are typically detected within 
rsoo- The rest are mostly point sources identified as back- 
ground AGNs. As it is difficult to distinguish AGNs from 
extended gas clumps within clusters at large radii (especially 
for ROSAT images), it is difficult to quantify the effects of 
clumps on the ellipticity measurements. At small radii, how- 
ever, detectable clumps are rare and do not change our results 
except in a few obvious cases. In order to robustly measure 
the global ICM shape while minimizing biases due to point 
sources and clumps, we mask them out and replace the holes 
with values from large-scale wavelet maps corresponding to 
scales greater than 32 pixels in each map. Improper treatment 
of point sources and clumps introduces significant differences 
in the ellipticity measurements. 
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The radial binning is set up logarithmically within rsoo such 
that the photon count in each bin is approximately the same. 
In addition, we ensure that there are at least 25 pixels in each 
radial bin. Except for the innermost bin in a few cluster, all 
bins have > 100 pixels. Finally, in our analysis, we allow the 
centroid to move for each annulus. Fixing the centroid will 
change the ellipticity typically by less than 10%. 

We estimate errors on elli pticity measureme nts using a 
bootstrap approach similar to Bu ote et ail (120021) . where we 
generate 100 realizations of the X-ray map and estimate the 
ellipticities using the method described above. The error is 
taken to be the dispersion of the bootstrap samples. The error 
is typically a few percent of the ellipticity value. 

We apply the same procedure above for our mock X-ray 
maps for consistent comparison with the observed clusters. 

3.3. Dynamical state 

To characterize the dynamical state of each cluster, we clas- 
sify our observational sample into relaxed and unrelaxed clus- 
ters. This can be accomplished in two ways. The first one 
is qualitative, where we classify clusters based upon visual 
examinations. If their X-ray morphologies are regular and 
have no secondary peaks or filaments, we classify them as 
relaxed, otherwise the y are unrelaxed (cf. Section 4.1.3 in 
Vikhli nin et al.ll2009al) . This is the classification that we used 
to classify our mock images of simulated clusters in Paper I. 

The second method is more quantitative: we classify a 
cluster as relaxed, if the logarithmic slope of its gas den- 
sity at r = 0.04rsno is les s than -0.5 for the Chandra sam- 
ple (IVikhlinin et al] 120071) . We were not able to apply this 
method to our CSF sample because their gas density slopes 
are too steep at r < 0.1 rsoo due to "overcooling." For a fair 
comparison between observation and simulations, we classify 
the cluster dynamical state visually. We have checked that the 
two methods result in similar samples of relaxed clusters, and 
our results for average ellipticity are within l-cr for the two 
observational samples defined in these ways. For the compar- 
ison with simulated cluster ellipticity profiles, we focus our 
comparison only on relaxed clusters because in this case el- 
lipticity primarily reflects ICM physics rather than spurious 
effects due to transient substructure or incomplete relaxation 
of the gas distribution during and after mergers. 

4. RESULTS 

4.1. Test cases: A2597, MKW 3s, and Hydra A 

To demonstrate the robustness of our method applied to ob- 
served clusters, we first examine the ellipticity measurements 
of three relaxed, cool-core clusters, A2597, MKW 3s, and 
Hydra A, for which indepe ndent cluster shap e measurements 
exist in the recent literature (Fang et al. 2009). Note, however, 
that all three clusters exhibit AGN activity in their central re- 
gions. 

As an illustration, we show in Figure [2] the X-ray image of 
A2597 observed with Chandra and ROSAT. On each image, 
we show the isophotal contours and the fitted ellipse from our 
ellipticity estimation. There is a good agreement between the 
fitted ellipses and the X-ray contours, demonstrating the ro- 
bustness of our method of estimating ellipticity on observed 
X-ray images of clusters. Note that there are X-ray cavities 
associated with AGN a ctivity in the central region of A2597 
dMcNamara et al.ll2001l) . However, structures associated with 
AGN activity are not discernible in most of our images be- 
cause of lack of resolution, and therefore our fitted ellipses 
are not affected by these structures. 



TABLE 2 

Mean Ellipticity of thez < 0.1 Chandra and ROSAT Clusters 



r/'"500 = 


0.05 


0.1 


0.3 


Relaxed 


0.17±0.03 


0.18 ±0.03 


0.20 ±0.02 


Unrelaxed 


0.38 ±0.20 


0.28 ±0.06 


0.23 ±0.09 


All 


0.19±0.06 


0.19±0.03 


0.21 ±0.03 



Figure [3] presents a comparison of the ellipticity profiles of 
A2597, MKW 3s, and Hydra A. The triangles and circles are 
the ellipticity measurement of our Chandra and ROSAT im- 
ages respectively. Crosses are independe nt analyses o f Cha n- 
dra and ROSAT/PSPC observations by iFang et alJ ((2009)0 
Here, we computed the ellipticities using all pixels within a 
given radius rather than pixels in annular shells, in order to 
compare to the results of Fan g et ail ([2009). Figure [3] shows 
that our ellipticity profile measurem ents are in good agree- 
ment with the shape measurements of lFang et alJ ?2009). 

4.2. Comparison of Mean Ellipticity Profiles in Local X-Ray 
Clusters vs. ACDM Clusters 

Figured] shows the average ellipticity profile for the 31 re- 
laxed clusters with both Chandra and ROSAT data (see Ta- 
ble Q]). The mean ellipticity of observed clusters (defined 
in Equation ( [3])) is e k 0.2 ±0.1 across the radial range 
0.05 < r/Vsoo < 1. The error bars show the jackknife errors 
on the error- weighted mean ellipticities. 

In the same figure we compare the observed mean ellip- 
ticity profile to the profiles of the relaxed NR and CSF clus- 
ters, shown in the left and right panels, respectively. Note 
that for the simulated clusters, we do not report ellipticity for 
r < 0.03r5oo because it reaches the limit of the spatial resolu- 
tion. The simulated ellipticity profiles are shown separately 
for each x,y,z projection. Note that the x,y,z axes correspond 
to the coordinates of the simulation box, and the projection is 
done along each orthogonal axis. Since the orientation of the 
cluster is arbitrary with respect to the box, each projection is 
a random sight line through each cluster. 

The number of relaxed clusters in our simulated sample is 
fairly small, so we also show the ellipticity profiles averaged 
over x, y and z projections independently (shown as solid, 
dashed, and dot-dashed lines, respectively) for both NR and 
CSF runs. This illustrates that the general trends seen in the 
NR and CSF are not significantly affected by orientation ef- 
fects due to small number statistics, as similar trends are seen 
in each independent projection. 

In the core regions (r < 0.1 rsoo), the NR ellipticities de- 
scribe the observations relatively well, with NR ellipticities 
close to £nr k 0.2. On the other hand, the CSF ellipticities 
are significantly higher than the observed ellipticities in the 
core regions, with ellipticities of ecsf ~ 0.3-0.5 due to the 
cold gas slab forming in the center via a cooling flow. 

At larger radii, 0.1 < r/r^a < 1, the ellipticity profiles of 
the NR clusters continue to match the observed ellipticity pro- 
files quite well with e w £nr ~ 0.2, although there is a hint 
that the simulated NR clusters may have slightly lower val- 
ues of e than the observed clusters. On the other hand, the 
ellipticity profiles in the CSF simulation clusters are consid- 
erably lower than the observed ellipticities over this range of 
radii. For example, at r = 0.2rsoo, £csf ~ 0.13, compared to 

8 Note that we have corrected for the radius used in Fang et al. 12009) by 
normalizing the radial scale using the same rjoo reported in TablefT] 
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FIG. 2. — X-ray image of A2597 observed by ROSAT (left panel) and Chandra (right panel. Shown in both panels are the X-ray isophotal contours (red thin 
lines) and the fitted ellipses (green dashed lines). The pixel size is 5.97 kpc for the Chandra image and 22.7 kpc for the ROSAT image. The size of the major axis 
of the largest ellipse is 0.065r5oo in the Chandra image and 0.28r5oo in the ROSAT image, where rsoo = 948 kpc for this cluster. Note that the brightness scales 
of the two images are different. 
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FIG. 3. — Comparison of ellipticity profiles of three relaxed, cool-core clusters: A2597, MKW 3s and Hyrda A. Here we compute the cumulative ellipticity using 
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Crosses indicate the ellipticity measurements from Fang et al. 12009), which also combine Chandra data at small radii and ROSAT/PSPC observations at larger 
radii. Error bars are estimated using 100 bootstrap samples. The size of the circles indicates the size of the largest error bar. 



e « 0.22 for the observed clusters. As discussed in the intro- 
duction, the differences between the ellipticities of the NR and 
CSF clusters at large radii reflect the more spherical shapes 
of the overall cluster potentials, and correspondingly more 
spherical shapes of the ICM in the CSF simulations, due to 
the formation of the massive central galaxies in those clus- 
ters. Observations thus clearly indicate that the central galaxy 
in the CSF simulations is too massive compared to observa- 
tions. This provides independent evidence for overcooling in 
the CSF simulations. 

To investigate the dependence of ellipticity on cluster 
mass, we show cluster ellipticities measured at r/rsoo = 
0.05,0.1,anc/0.3 for the observed clusters as a function of 
M500 in Figure The average ellipticity at r = 0.05rsoo is 
derived from the Chandra observations, while the ellipticity 
values at the larger radii are derived from the ROSAT sample. 



There is no obvious dependence of ellipticity on cluster mass, 
consistent with the simulation results presented in our Paper I. 
The scatter in ellipticity among clusters is also similar across 
all radii considered. A weak trend of ellipticity with cluster 
mass may be present, but this sample of clusters is too small 
to constrain it robustly. 

4.3. Effects of Dissipation on Ellipticity Profile 

To investigate the dependence of the ellipticity on baryon 
dissipation, we re-simulated one of the clusters in simula- 
tion sample, CL7, with gas cooling, star formation, and stellar 
feedback turned off atz < 2.0 (CSFoff). CL7 is an isolated, re- 
laxed cluster at z = 0, which did not experience a major merger 
for several Gyr prior to z = 0. The CSFoff experiment resulted 
in a final cluster with a stellar mass within r = 50 kpc about 2 
times smaller than the stellar mass in the corresponding CSF 
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simulation. 

We estimate the ellipticity for CL7 using the same methods 
that we applied to the observed clusters. Figure [6] shows the 
ellipticity profiles for the three orthogonal projections of CL7 
(left, middle, and right panels for the x,y, and z projections, 
respectively) for the three runs. The red (dot-dashed), green 
(dashed) and blue (solid) lines represent the ellipticity profiles 
for the NR, CSFoff and CSF runs respectively. The agree- 
ment between the CSF and CSFoff profiles for r > O.lrsoo 
shows that the effect of baryon condensation on the shape of 
the potential and ICM does not scale linearly with stellar mass 
within rsoo- Evaluation of this dependence will require further 
investigation with a much larger sample of clusters. 

Note that for r < 0.1 rsoo there is a significant difference 
between CSF and CSFoff runs. The absence of cooling after 
z = 2 prevented the formation of a cooling flow and extended 



flattened slab in the central gas distribution. The absence of 
this gaseous slab resulted in ellipticity values at r < 0.1 rsoo 
closer to the values in the outer regions of the cluster. This 
clearly illustrates that the oblate slab in the core of this cluster 
in the CSF run is due to cooling during later stages of cluster 
evolution. 

5. SUMMARY AND DISCUSSION 

In this study, we have compared mean ellipticity profiles 
of a sample of clusters simulated in the concordance ACDM 
cosmology with ellipticity profiles measured with Chandra 
and ROSAT X-ray observations of galaxy clusters. 

Our comparison presents several improvements over the 
previous work. In terms of analysis method, we have used 
mock X-ray photon maps of simulated clusters and used the 
same analysis pipeline on both the mock and real maps. This 
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FIG. 6. — The ellipticity profiles for the three orthogonal projections of CL7 (left, middle, and right panels for the x,y, and z projections, respectively) for 
the three runs with different physics. The red (dot-dashed), blue (solid), and green (dashed) lines represent the ellipticity profiles for the NR, CSF, and CSFoff 
(cooling and star formation turned off at redshift z < 2) runs respectively. 



enabled more direct and consistent comparison between simu- 
lations and observations. In terms of sample size, we have sig- 
nificantly expanded the size of observat ional sample by ana- 
ly zing 3 1 local clus ters compared to 2 dFlores et alj2007l) and 
9 dFang et al. 2009). Using this sample, we measured the av- 
erage ellipticity of local relaxed clusters to be e = 0.18 ± 0.05 
in the radial range of 0.04 < r/r 5 oo < 1. Our results indicate 
that the shapes of X-ray emitting gas in galaxy clusters, espe- 
cially at large radii, can be used to place constraints on cluster 
gas physics, and are potential probes of the history of bary- 
onic cooling in galaxy clusters. Key results of the paper are 
summarized below. 

We showed that clusters formed in simulations with gas 
cooling and star formation (the CSF runs) generally predict 
ellipticities considerably smaller (i.e., a rounder ICM) than 
the observed clusters at radii 0.1 < r/rsoo < 1, but they pre- 
dict much higher ellipticities than the observed clusters in the 
core region (r < QAr^oo). The non-radiative simulations, on 
the other hand, have ellipticity profiles much closer to the ob- 
served cluster ellipticities in the range r > 0.05r5oo. 

The amount of baryonic dissipation occurring during clus- 
ter formation controls the stellar masses of central cluster 
galaxies. The stellar masses of central cluster galaxies, in 
turn, control the shapes of the global potentials and the 
ICM of clusters at large cluster-centric radii. Our compari- 
son therefore suggests that the masses of central galaxies in 
observed clusters are much smaller than the corresponding 
masses in our CSF simulations with gas cooling and star for- 
mation. This is a manifestation and an independent confirma- 
tion of the "overcooling" commonly seen in cluster formation 
simulations with radiative cooling and star formation (e.g., 
Borgani & Kravtsov 2009). To probe the degree to which cen- 
tral galaxy mass is overestimated we have re-simulated one of 
our relaxed clusters with cooling and star formation turned off 
at z < 2. The cluster in this re-simulation had a considerably 
smaller final stellar mass, but still exhibited an ellipticity pro- 
file consistent with the CSF simulations at large radii. This 
indicates that the sensitivity of the ICM ellipticity to the mass 
of central galaxy is nonlinear and will need to be quantified 
using larger samples of objects. 

In the central regions of clusters, r < O.lrsoo, the ICM 
shapes in clusters formed in simulations with cooling are gen- 
erally oblate due to cooling flows that develop in the central 
regions and lead to the formation of rotating thick gas slabs. 



This result is consistent with the results of iFang et al.l (120091) . 
who analyzed a subset of the cluster sample we study in this 
paper. However, these authors found the sizes of the oblate 
slabs to be a bout a factor of t wo larger than we find here (see 
discussion in lLau etafll201 ll) . The formation of such cooling 
flows and slabs occurs during late stages of cluster evolution. 
For instance, in the CSFoff re-simulation, with cooling and 
star formation turned off at z < 2, the shape of the gas dis- 
tribution near the cluster center is not oblate and the cluster 
ellipticity is nearly constant with radius. 

The ICM shapes in the central and outer regions of clusters 
therefore probe and constrain different regimes of baryonic 
cooling. The shapes at large radii are sensitive to the over- 
all amount of gas that has cooled and condensed to form a 
central cluster galaxy throughout the entire cluster evolution. 
We reiterate that our results indicate that the relationship be- 
tween ICM shapes at large radii and the amount of gas that 
has cooled to form the central cluster galaxy is nonlinear and 
may be more well calibrated with a larger simulation set. Al- 
ternatively, ICM ellipticities in the central regions of clusters 
probe recent and ongoing gas cooling. 

As discussed in Paper I, random, isotropic gas motions in 
the cluster cores result in more spherical gas shapes in the cen- 
tral regions of non-radiative clusters. This leads to a downturn 
in the ellipticity profiles toward more spherical ellipticity val- 
ues at r/rsoo < 0.1 (see the left panel of Figure|4j. The mean 
ellipticity profiles of observed clusters do not exhibit such a 
downturn at these radii. This implies that real clusters do not 
have gas motions as significant as those in our NR simula- 
tions. The likely reason for this difference is that observed, re- 
laxed clustersjiavejrionotonically decreasing entropy profiles 
(e.g JCavagnolo et al. 2009), while the simulated NR clusters 
have flat entropy profiles at 0.15 < r/rsm < 0.2 which are 
more conducive to convective random gas motions. 

A recent work by lBiffi et al.l (1201 ll) also studied the effects 
of baryon cooling on gas ellipticity and velocity structure us- 
ing GADGET SPH code. While our results ba sed on the 
AMR grid code are consistent with the results of Bi ffi et all 
(2011) in the outer regions of clusters, there are noticeable 
differences in the central regions. Specifically, the cooling 
runs performed by the SPH code do not produce a rotation- 
ally supported gas disk in the central regions of relaxed clus- 
ters, as seen in the AMR simulations. One possible explana- 
tion for such disparate results is the difference in disruption of 
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gas c lumps in SPH and AMR simulations (e.g., Age rtz etldl 
2007). Gas clumps in the SPH code do not disrupt as effi- 
ciently as those in the AMR code due to SPH's poor ability 
to capture fluid instabilities (Kelvin-Helmholtz in particular). 
Since the unstripped gas loses angular momentum while the 
stripped gas does not, rotationally supported gas disk does not 
build up in the core of SPH clusters as extended as that in the 
AMR clusters. 

The ellipticity profiles of clusters in non-radiative simula- 
tions agree with obs erved ellip ti cities well, which is consis- 
tent with results of Fan g et all (120091) . who used a smaller 
sample of clusters with Chandra and ROSAT/PSPC data to 
de rive mean ell i pticity profile. Similar results were also found 
bv lFlores et al.l (120071) . who used a subset of our non-radiative 
cluster simulation and showed that their shapes are in general 
agreement with o bservations . Our r esults are also consistent 
with the results of Kawahara (2010), who showed that distri- 
bution of cluster ellipticities derived from XMM-Newton ob- 
servations are consistent with the predi ctions of diss i pation - 
less A^-body simulations of clusters by iJing & Sutol (120021) . 
The broad accord between our results and these previous stud- 
ies suggests that the ellipticity profiles in clusters formed in 
non-radiative simulations are in reasonable agreement with 
the ellipticity profiles of observed clusters. 

One possible concern with the comparison between simu- 
lations and observations presented in this paper is the pos- 
sible bias in the ellipticities of simulated clusters due to the 
fact that our cluster sample was simulated assuming erg = 0.9, 
while the most recent estimates indicate e ra = 0.80 ±0.02 (e.g ., 
IVikhhnin et aTll2009bt iRozo et al.ll20T51 Uarosiket al.ll2011l) . 
The effect of a lo wer a% is that halos form later, leading to 
higher ellipticities (Allg ood" et al.l l2~006). This effect is likely 
to be small, so that it would not change o ur qualitative conclu- 
sions. For example, Mac cio et al.1 (120081) show that effect of 
changing a% from 0.9 to 0.8 changes average minor-to-major 
axis ratios, c/a, of halos by only rs 0.03. This systematic 
shift is considerably smaller than the differences discussed in 
the present paper. Furthermore, the sense of this effect would 
be to bring the NR simulations into better accord with the ob- 
served cluster sample. 

On the observational side, a possible concern is that buoy- 
ant bubbles or jets associated with the central AGN activity 
could modify the morphology of gas in the core regions, re- 
sulting in higher observed ellipticities. We have visually ex- 



amined our observed sample for any signatures of AGN ac- 
tivity and have excluded regions that seem to be affected by 
AGNs. However, it is still possible that some regions influ- 
enced by AGN activity are unresolved in our current sample, 
or that buoyant bubbles have had an impact on ICM morphol- 
ogy, but subsequently subsided so that they are not observed 
at present. 

In this paper we compared the shapes of simulated ACDM 
clusters to the ICM observed in a local X-ray cluster sample. 
A natural extension of this work is to probe the ICM shape of 
higher redshift clusters with deeper X-ray imaging and high- 
resolution Sunyaev-Zel'dovich (SZ) observations. It would 
also be useful to carry out such comparisons with larger sam- 
ples of simulated clusters for a range of cooling and feedback 
models, resulting in a wide and well-sampled range of central 
galaxy stellar masses. High-resolution SZ imaging is a poten- 
tially powerful means of extending this analysis to higher red- 
shifts, because of the redshift independence of the SZ effect. 
Several high-resolution SZ imaging experiments are in oper- 
ation or under construction, including ALMA0, CARMAEH 
CCAlEi and MUSTANCQ Detailed comparisons between 
SZ shapes in simulations and observations should provide fur- 
ther tests of the baryon dissipation in galaxy clusters over a 
wide range of redshifts. Such comparisons should inform the- 
ories of the formation of galaxies and galaxy clusters and may 
play a complementary role in efforts to refine constraints on 
the concordance cosmological model. 
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